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We present improved Coulomb correction formulae for Bose-Einstein correlations including 
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also exchange term and use them to calculate appropriate correction factors for several source 
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functions. It is found that Coulomb correction to the exchange function in the Bose-Einstein 



correlations cannot be neglected. 
SULDP-1995-3 TU477 SINS- 1995-1 Introduction: Recently many experimental groups 

have investigated Bose-Einstein correlations (BEC) in high energy hadronic collisions [jj], e + e~ 
annihilations and heavy- ion collisions || 0] . The high quality of data obtained already (and ex- 



On 

a,: 
a: 

> 

X. 

pected in near future, especially for heavy-ion collisions) make the analysis of experimental results 
sensitive to all possible corrections, especially to those due to final state interactions because of 
the Coulomb interactions and the strong interactions |j| . In the present letter we shall concentrate 
on the Coulomb corrections only. 

Several authors have calculated Coulomb corrections to the BEC by their own methods |(| |7[ 
|§|, ||, |l(| but this problem is still controversial jn], [ll[] . In this letter we shall concentrate on the 
formula provided by Bowler fl0|| , because it includes all orders of the parameter r\ (defined below). 
He demonstrated that widely used zero range Gamow factor substantially overestimates Coulomb 
corrections (and hence also the true magnitude of BEC). However, Bowler did not calculate the 
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Coulomb correction to the exchange function present in BEC formulae. We investigate therefore 
this problem and in particular: (1) we calculate the Coulomb corrections using the exact formula 
(to all orders in parameter 77) and (2) we use its approximate form (retaining only terms linear 
in rf) to calculate explicit Coulomb wave function in this approximation |?], |l2], [l3| and with its 
help we provide formulae for Coulomb corrections for the exchange function in the BEC for several 



source functions and compare some of them with those presented in ref . 1 10 . We demonstrate that 
the presence of exchange term diminishes correction factor even more than anticipated in |l0| ] and 
therefore it should be included in analysis of experimental data. 

Theoretical calculation of BEC with Coulomb wave function: To write down an amplitude 
A12 satisfying Bose-Einstein statistics it is convenient to decompose the wave function of identical 
(charged in our case) bosons with momenta p\ and P2 into the wave function of the center-of-mass 
system (cm.) with total momentum P = \(p\ + P2) an d the inner wave function with relative 
momentum Q — (pi~p2) = 2q. It allows us to express A12 in terms of the confluent hypergeometric 
function $ Q: 

A U = ^[*(q,r) + * s (q,r)], (1) 
#(q,r) = r(l + ir))e- vv/2 e iclT $(-iri;l;iqr(l-cos6)) , 
*s(q,r) = r(l+^)e- 7r? ' /2 e- iq - r $(-i7 ? ;l;i g r(l + cos6')), 

where r — x\ — X2 and the parameter 7/ = ma/2q. Using now the Kummer's first formula for the 
confluent hypergeometric functions that appears in the cross term of | A12 | , 

$(a; 7; z) = e z $(7 - a; 7; -z), 

we calculate first the exact formula of Coulomb correction (i.e., the one that is exact to all orders 
of parameter rf) including also the exchange function in BEC. Namely, assuming factorization in 
the source functions, p{r{)p{r2) = p(R)p(r) (here R = \{x\ + X2)), one obtains the following 
expression for theoretical BEC formula: 

N ±±/ N BG = J_J p{R)d 3 R J p ( r)d 3 r | Ai2 |2 
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(n + m)\ \ irj J \ irj 
where G(q) — 2irr]/(e 2mi — 1) is the Gamow factor and 



''.■<"•'"> = ^ / drr 2 +^p(r), A n = + n) -L 
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The second term in the squared parenthesis in eq.(J2|) is a new term which is due to the exchange 
function in the BEC. The above expression represents our main result. However, it has disadvantage 
that it is difficult to separate in it the exchange function and its Coulomb correction. Therefore, 
in order to know the separate contributions we must either subtract from it the exchange function 
or to use the approximate formula for the Coulomb correction. To calculate it we expand function 
<3? in powers of r\ and retaining only linear terms we get 

$(-mj; l;ix) = 1 + rjSi(x) - irj(Ci(x) — 7 B — ln(a:)) + 0(r? 2 ) , 
$(-Mj; l;ix) — 1 + r]Si(x) - ir](Ci(x) - 7 E - ln(x)) + 0(r? 2 ) 

where x — iqr(l — cos8), x = iqr(l — cos9), "f E is the Euler's number and Si(x) and Ci(ir)— "f E — ln(x) 
are the sine and cosine integral, respectively, defined as: 

^ (2n + l)!(2n + l) ^ (2n)!(2n) 

With such approximated $'s we obtain the following approximate formula for BEC with Coulomb 
corrections included: 



n (±±) /n bg = J_J p{R)d 3 R J p ( r )d 3 r\A 12 \ 



2 



where (A = 2qr) 
h(q) = 



G(q) 

= h(q)+h{q) (3) 
= 1 + Sic + Sec + E2B 

= (1+< ' c+fe) ( 1 + Trra?) < 4 » 



4. / P{r y d r { 1 + (2 ( w+ ) i ;y n + 1) /_ x (l - cos.)— .cos. 



1 + 4tt • 2t? / p(r)r 2 dr V 

J ^ (2n+l 



(_l)M 2 ™ +1 



^ (2n+ l)!(2n+ l)(2n + 2) 



1 + *io , (5) 
47r J p{r)r 2 dr |~~^ l~ ^[SpGzr) + Cp(gr)]l 

S 2B + S EC , (6) 



and 



~ (-l) 2 «+ 1 (gr) 2 "+ 1 6(2n+l) 
Sp(gr) = 2^ - 



n=0 



(2n + l)!(2n+ 1) 



with 

0(2n+l) = J {1 - cos 9) 2n+1 cos(A cos 9)d cos 9 , 

Q(2n) = T J {1± cos 9) 2n sin(A cos 9)d cos 9 . (8) 

(Notice that 6(1) = 2z*A, 0(3) = ^-^+^^—^^,...axidQ(2) = ^-^A-^A, . ..). 
The normalization-like factor (1 + Sic + Sec) m e 1-(@) differs from that in ref. only by the 
additional exchange term Sec- There is therefore the following (approximate) relation between 
our correction factor, C° urs , and that of Bowler pTOj ] , C l 



~i Bowler . 



c Ours = (jBowler + ^ . (g) 

where G{q) is Gamow factor. 

Coulomb corrections and source functions: In order to provide numerical estimations of our 
new correction factor we calculate analytical expressions for necessary ingredients of both exact 
result as given by eq.(|J), Ir, and the approximate one provided by eqs.([|) and (^), Sf c , E 2 b and 
Sec f° r several different source functions (superscripts below denote the type of source considered): 
(SI) Exponential source function, p(r) = |^ exp(— fir): 

Taking into account only leading terms in 0(2n + 1) and 0(2n) in eqs. (pj) one has: 

/|( „, m) _ (i)" + '"<H±^, oo, 

n=0 

E ' B = (1 + J2q/W? ' (12) 



+ 1 f ( ~ i r„ (2 ; 1 + 2) (i + g/OT-^ ^ + 3) Mctan w)) • (i3 » 

Our results for this source function are given in Fig. 1. Notice that there is significant systematic 
difference between C Bowler and C° urs correction terms whereas C° urs is practically the same 
whenever calculated using exact formula (||) or approximate one (^) (in the range of q considered 
the differences are of the order of 1%). 

(S2) Gaussian source distribution, p(r) = exp(— f3 2 r 2 ): 
For this type of source function we have: 

T G > \ 2 / 1 \ , _ N ^ / n + m + 1 
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2£(n,m) = ^[g] (n + m+l)T( ), (14) 



G 4 V » (-!)"(» + !)! /2 g \ 2 " +1 

1C " ^^ (2n + 2)!(2n + l) UJ ' ( } 

= exp f-^J , (16) 

2,, / ,'\f (-l)»(n+l)l) /2<A 2 " +1 / 1 3 ? 2 \ 

In Fig. 2 we show typical examples of our correction term for Gaussian source distribution. Also 
here exact and approximate formulas lead practically to the same results at lower values of q 
(compatible with those in Fig. 1) but differ substantially for larger q's. The difference between 
our results and the result obtained using Bowler's formula (cf. eq.(|27]) below) exists also here. 
(S3) Modified Bessel source functions, p(r) — J^KoiPr) and p(r) = -^rKi^r): 
In the case of source functions described by the Modified Bessel functions Kq or K\ |H| [l(| we 
have the following expressions: 
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I^n,m) = —[j) r^^jr^^j, (18) 
s k 8,^ (-!)»{(» + I)!} 2 /4 g \ 2 " +1 

^ " 7r ^ o (2n + 2)!(2n + l)UJ ' ( j 

E ™ = (i + WWW ' (20) 
s k ±i v (-!)" f^] 2 ^ 1 v ( ' 1)m{(n + m)!}2 % " 

2r? ^ (-1)" /4 g \ 2 " +1 ^ (-ir{(n + rn+l)!} 2 /2g\ 2 "' 

V5F^(2n+l)!(2n + l) U>/ ^ r(m+|)m! ^ J ' [ 1 

and 
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The correction factors for these source functions are presented in Fig. 3 which closely resembles 
Fig. la (again, approximate formula practically does not differ from the exact one whereas both 
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differ substantially from that calculated according to [[To|). 

Concluding remarks: We have calculated the exact and approximate formulae for Coulomb 
corrections used in the BEC (including exchange term) and specified them to the Exponential, 
Gaussian and Modified Bessel source functions. The difference between exact and approximate 
formulas (eq. (|^) and eq.(^), respectively) turns out to be smaller than 1 % (at least for q < 0.2 
GeV/c). However, due to the presence of the new correction for the exchange term (Sec above) 
our results are systematicall lower than those derived in pp| . To visualize it better we check if 
relation 

= 1 (26) 

OEC 

holds; E 2B and Sic above are quantities used in ref. [fTo| : 

iV±±/7V BG = (1 + Sic)(l + E 2B ). (27) 

It is clear from Fig. 4 that (|2^) does not hold (here for the standard value of /3 = 0.2 GeV used 
also in JIo| but we have checked that the same is true for the whole possible range of this param- 
eter). We conclude that correction term Sec m cc l s -(||) an d (0) cannot be neglected and that in 
calculations of Coulomb corrections eq.(p7|) should be replaced by eq.(||) (or eq.(||)). 

In ref. p0| the relative changes of the normalised Fourier transforms of a given source function, 



, were also introduced and investigated as yet another estimate of the importance of the Coulomb 
corrections: 

Ap _ (1 + 5 1C )(1 + E 2B ) - 1 

T - (l + E 2B )-l 1; (28) 

for the ideal BEC and 

A/3 fi x(l + S w )(l + E 2B ) + G-Hl-x)-l _ 

p x[(l+E 2B )-l] (y> 

for BEC containing contributions of the long lived resonances {L} — {77,^,77 ; c, b} with 1 — x 
denoting the fraction of pairs involving a daughter of {L}. In Figs. 5 and 6, using the exact formulae 
containing also the correction term Sec, we compare our results for ^ and ^ calculated for 
Exponential and Gaussian sources for different /3's (in Fig. 5) and different x's (in Fig. 6). Contrary 
to the previous cases, for these quantities the introduction of Coulomb correction to the exchange 
term in BEC practically does not change the results obtained by using eq.(|27]) from [0, due to 
the large denominator. 
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Figure Captions 

Fig. 1. Comparison of C Bowler with C rs for exponential source function and for different 
choices of parameter 0: (a) (3 = 0.2 GeV. (b) (3 = 0.1 GeV. 

Fig. 2. The same as in Fig. 1 but for Gaussian source distribution. 

Fig. 3. The same as in Fig. 1 but for Modified Bessel source functions (Kq for (a) and K\ for 
(b)) - only for j3 — 0.2 GeV in both cases. 

Fig. 4. Examinations of eq. (|2q) for different choices of source functions. 



Fig. 5. Aj5/p (cf. eq.(p8|)) for different f3's and for Exponential and Gaussian source distributions. 
In our calculations the exact formulae are used. 



Fig. 6. Ap R /p (cf. cq.(29)) for Exponential and Gaussian source functions and for different 
values of the fraction of long lived resonances parameter x. In our calculations the exact 
formulae are used. 
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